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Abstract 

Chaotic diffusion on periodic orbits (POs) is studied for the perturbed Arnol'd 
cat map on a cylinder, in a range of perturbation parameters corresponding 
to an extended structural-stability regime of the system on the torus. The 
diffusion coefficient is calculated using the following PO formulas: (a) The 
curvature expansion of the Ruelle zeta function, (b) The average of the PO 
winding-number squared, w 2 , weighted by a stability factor, (c) The uni- 
form (nonweighted) average of w 2 . The results from formulas (a) and (b) 
agree very well with those obtained by standard methods, for all the pertur- 
bation parameters considered. Formula (c) gives reasonably accurate results 
for sufficiently small parameters corresponding also to cases of a considerably 
nonuniform hyperbolicity. This is due to uniformity sum rules satisfied by 
the PO Lyapunov eigenvalues at fixed w. These sum rules follow from general 
arguments and are supported by much numerical evidence. 

PACS numbers: 05.45.Ac, 05.45.Mt, 45.05. +x 
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I. INTRODUCTION 



Understanding to what extent chaotic motion in Hamiltonian systems exhibits random 
properties such as diffusion is a problem of both fundamental and practical importance. 
The existence of deterministic chaotic diffusion has been approximately established using a 
variety of approaches [1-20]. A systematic approach is based on the hierarchy of periodic 
orbits (POs) embedded in the chaotic region [10-20]. Let us summarize the main ideas of 
this approach by considering, for definiteness, the kicked-rotor maps on the cylinder 

(1) 

where / is angular momentum, x is angle, and the force function f(x) satisfies f(x+l) = f(x) 
and f(—x) = —f(x). The diffusion coefficient for (1) is formally defined by: 

D = lim D £ (m), D s {m) = «^ ~ ^ ( 2 ) 

m— >oo ZlTL 

where ( )s denotes average over an ensemble £ = {(xo, /o)} of initial conditions in a chaotic 
region. Well known problems in a reliable numerical calculation of D are: (a) The roundoff 
errors caused by the chaotic exponential instability, (b) The ambiguity in the choice of £ 
and the iteration time m. These problems are systematically solved [10] by choosing £ as 
the ensemble U n of all the primitive POs of period n in the chaotic region and m = n. For 
the map (1), a PO p G U n is generally defined by initial conditions (x^\ l^) satisfying 

lip) _ lip) I 7 „ T ip) _ JP) (o\ 

where n is the smallest integer for which (3) holds with integer w p , the winding number of 
PO p. Since (1) is essentially periodic in (x, I) with period 1, {x^\ l^) can be restricted 
to a unit torus, —0.5 < x^\ < 0.5. The ensemble U n may be viewed as an invariant 
"level n" approximation of the chaotic region and the diffusion rate on U n is given by [10]: 
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where N(n) is the total number of POs in U n and N w {n) is the number of POs in U n with 
w p = w. Similarly, one can associate diffusion rates with subensembles of U n having well- 
defined dynamical characteristics [10,11]. For uniformly hyperbolic systems, (4) is expected 
to approximate well the diffusion coefficient (2) associated with generic ensembles of aperi- 
odic chaotic orbits. In fact, in the case of the cat and sawtooth maps [2], (4) gives the exact 
value of D for the cat maps and approximates very well D for the sawtooth maps [10]. 



The diffusion coefficient for generic chaotic ensembles in hyperbolic systems is given by 
the exact PO formula [16,19,21]: 

1 s)/d(3 2 



D 



2 d(-i((3, s) /da 
Here s) is the Ruelle zeta function [21], 



(5) 

(3=s=0 



C~\P, s) = l[[l-exp(f3w p -sn p )\A p r l ], (6) 

p 

where the product is over all the primitive POs, n p is the period of PO p, and A p is the 
associated Lyapunov eigenvalue (|A P | > 1). One can express (6) as a power series in exp(— s): 

oo 

C\f3, s) = l + Y^c n (P)exp(-sn), (7) 

n=l 

where the n > 1 terms are known as "curvatures" [19]. The convergence of (7) is generally 
better than that of (6). Formula (5) with a (truncated) curvature expansion (7) has been 
applied to several systems [13,14,16-20]. For uniformly hyperbolic maps with a complete 
symbolic dynamics (i.e., whose grammar is unrestricted by "pruning rules"), all the curva- 
ture terms in (7) vanish identically and the exact value of D from (5) coincides precisely 
with the diffusion rate (4) for n — 1. Examples of such trivial systems are ID piecewise 
linear maps [13,14] and chains of coupled baker maps [16,21]. The application of (5) and 
other PO formulas [15] to a more realistic system, the periodic Lorentz gas, gives results [18] 
that are within 8% of the values of D obtained by standard methods. For standard maps 
(1), the quasilinear (strong-chaos) limit of D is approximated by using just POs of period 



n = 1 and n = 2 [17]. For the cat and sawtooth maps, the quasilinear approximation of D 
is reproduced by formulas related to (5) [20]. 

In this paper, chaotic diffusion on POs is studied for a nontrivial Hamiltonian system 
exhibiting a transition from uniform to nonuniform hyperbolicity as a parameter is varied. 
This is the perturbed Arnol'd cat map on the cylinder, defined by (1) with 

f(x) = f (x) + ^sm(2nx), (8) 

where fo(x) = x for \x\ < 0.5, / ( — 0.5) = 0, f (x + 1) = fo(x), and k is a perturbation 
parameter. This system, with the definition / (— 0.5) = —0.5, is usually considered on a 
torus, —0.5 < x, I < 0.5. Perturbed cat maps on the torus have attracted much attention 
recently in the context of "quantum chaos" [22-24]. Anosov theorem [25] states that the 
dynamics on the torus for sufficiently small k, k < k c , is topologically equivalent to that of 
the unperturbed (k = 0) system (in particular, the system is completely chaotic for k < k c ). 
This expresses the well known structural stability of cat maps (see more details in Sec. II, 
where we determine k c ~ 0.437). Actually, we provide numerical evidence in Sec. II that 
the structural-stability regime extends, at least approximately, beyond k c , up to k ~ 1. A 
fully chaotic regime is observed up to k ~ 1.5. For larger perturbations, stability islands 
born by bifurcation, leading to a significant mixed phase space for k > 1.7 (see Fig. 1). We 
emphasize that unperturbed cat maps already feature a very nontrivial symbolic dynamics 
with nonexplicit pruning rules given by an infinite set of inequalities [26]. As a result, all the 
curvature terms in zeta-function expansions are nonvanishing. See, e.g., an exact expression 
of C _1 (0, s) for cat maps derived in Sec. II. The relevant dynamics on the cylinder can be 
easily inferred from that on the torus. 

In Sec. Ill, we calculate accurately the diffusion coefficient D for k up to k ps 1 using 
several PO formulas: (a) The curvature-expansion formula [(5) with (7)]. (b) The average 
of Wp weighted by a stability factor. This formula was used in Refs. [15,18] and we give 
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a derivation of it from (5). (c) The nonweighted-average formula (4). The convergence 
of each formula as the order n of approximation is increased is verified. When compared 
with results obtained by standard methods, the results from formulas (a) and (b) appear 
to be more accurate than those in previous works. The relative difference between the PO 
and standard results is not larger than 0.4% when formula (b) is used and not larger than 
1.7% when formula (a) is used. Formula (c) gives reasonably accurate results (within 2% 
of the standard results) for sufficiently small values of k corresponding also to cases of a 
considerably nonuniform hyperbolicity. This is due to uniformity sum rules satisfied by the 
PO Lyapunov eigenvalues at fixed winding number w. These sum rules follow from general 
arguments and are supported by much numerical evidence. 



II. CAT MAPS AND STRUCTURAL STABILITY 

Consider the unit torus T 2 : —0.5 < x, I < 0.5 and let z = (x, I). Hyperbolic cat maps 
on T 2 are defined by the map 0o : z' = A ■ z mod T 2 , where A is a 2 x 2 integer matrix 
with det(A) = 1 and Tr(A) > 2. While these maps are uniformly hyperbolic, they feature a 
very nontrivial symbolic dynamics [26]. As a result, all the curvature terms in zeta-function 
expansions are nonvanishing. We show this here by deriving an exact expression of C _1 (0, s) 
for O . First, the uniform hyperbolicity implies that the Lyapunov eigenvalue A p of any PO 
with period n is given by A p = A™, where A is the largest eigenvalue of A. Then, from (6), 

oo 

C\0, s) = U{ 1 -P n ^ N(n) > ( 9 ) 

n=l 

where p(s) = A _1 exp(— s) and N(n) is the number of primitive POs of period n. For 
p(s) < 1, we find from Rel. (9) that 

oo 

ln [C _1 (0, s)} = Y,N(j)^[l-p\s)} = 

-f>wE^ = -E^5>w). no) 

j=l i=l J n=l j\ n 
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where j \ n means that the positive integer j divides n. We now use the general relation 

J2jN(j) = p(u), (11) 

where P{n) is the number of periodic points of period n. For the cat maps one has [28] 

P(n) = \Tr{A n ) -2| = A n + A~ n -2. (12) 



Using (11) and (12) in (10), we obtain, for s > 0, 

[1 — exp(— s)] [1 — A~ 2 exp(— s) 

r yii) = 

n 

n=l 



C\0, s) = exp 



[1- A- 1 exp(-s)] 2 



(13) 



The explicit curvature expansion of ( x (0, s) for s > is easily found from (13): 

oo 

(~\0, s) = 1 -P(l)^nA- n exp(-sn). (14) 

n=l 

Thus, all the curvatures c n (0), n > 1, are nonvanishing. 

A perturbed cat map on T 2 is given by 4> K : z' = A ■ z + kF(z) mod T 2 , where 
F(z) is a smooth vector field periodic on T 2 . Anosov theorem [25] states that for suffi- 
ciently small k,, k < k c , <p K is topologically conjugate to O by a continuous near-identity 
map H K , <p K = H K o O o H^ 1 . Thus, any orbit O k of (f> K , in particular a PO, can 
be written as O k = H K Oo, where Oo is some orbit of 0o- We emphasize that the re- 
lation (p K = H K o (f) o H~ x does not imply that the Lyapunov eigenvalue A p of a PO 
O k = H K O is equal to that of O since the map H K is not differentiate [25]. The bound 
k c is generally determined by the inequality: max zgT 2 (|/t<9F/<9z • z|/|z|) < 1 — A -1 , where 
|z| = (x 2 + I 2 ) 1 / 2 . Choosing A = (2, 1; 1, 1) (corresponding to the Arnol'd cat map [25]) 
and F(z) = (l/27r) sin(27ra;)(l, 1), we see that <p K is just the map (1) with (8) and the 
definition / ( — 0.5) = —0.5 [instead of f (— 0.5) = 0]. In this case, A = (3 + \/5)/2 and 
\d¥/dz ■ z| = V2\x\ cos(27ra;). From the inequality above we then get k < k c m 0.437. 



The POs of a perturbed cat map for k < k c can be calculated very accurately as fol- 
lows. First, the POs O of the cat map are determined exactly using the techniques in Ref. 
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[27]. The perturbed POs O k = H K O are then computed by applying to O the map H K 
constructed iteratively from the nonlinear functional equation satisfied by it, starting from 
the solution of the "homological equation" (see details in the proof of Anosov theorem 
[25]). In this way, we have calculated all the POs O k = H K O of the perturbed Arnol'd cat 
map with periods n < 14 with an accuracy of at least 1CT 10 ; this accuracy was checked by 
direct iteration of the map. As a matter of fact, we found no problems in extending these 
calculations beyond k c ks 0.437. We have thus verified that all the POs K = H K O with 
n < 14 exist at least up to k — 0.946. In order to check to what extent these POs cover the 
set of all POs with n < 14 for k c < k < 0.946, we have used the C _1 (0, 0) = test [16-19]: 
the curvature expansion (7) for C _1 (0, 0), restricted to the set of POs O k = H K O , was 
calculated up to order n — 1, 14 for six values of k uniformly distributed in the interval 
< k < 0.946. The results are shown in Table I. We see that the behavior of C _1 (0, 0) as 
n increases is basically the same for both k > k c and k < k c , and there is a general trend 
of |C — 1 (0, 0)| to decrease. This indicates that the POs O k = H K O form, at least, a large 
fraction of all the POs for k c < k < 0.946, i.e., almost no bifurcations take place in this 
interval. Thus, the structural-stability regime appears to extend, at least approximately, up 
to k = 0.946. On the basis of this observation, only the POs O k = H K O were used in our 
calculations for < k < 0.946. 

On the cylinder, the system is described by the map (1) with (8) and / ( — 0.5) = 0. 
The discontinuity of (8) at x = —0.5 can be viewed as an infinite-slope (vertical) segment. 
Is is then easy to show that POs with points on x = —0.5 must have an infinite value of 
A p . Thus, they will not contribute to (6) and will not be considered. Clearly, the initial 
conditions for the relevant POs (having no point on x = —0.5) can be chosen as points of 
torus POs O k = H K O lying in the domain \x\, |/| < 0.5. The number N(n) of relevant POs 
used in our diffusion calculations in Sec. Ill is listed in Table II. The winding number w p is 
calculated from (3). Strictly speaking, there is no structural stability on the cylinder since 
w p generally changes when k is varied in any interval, in particular (0, k c ). 
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III. RESULTS FOR D AND UNIFORMITY SUM RULES 



We start by deriving a general formula [Eq. (17) below] giving an approximation to D in 
terms of all the primitive POs of period n. Let us calculate the derivatives <9 2 C _1 (/3, s)/d/3 2 
and <9C _1 (/5, s)/ds in (5) for (5 = and s > by direct differentiation of the infinite product 
(6), which is convergent and nonvanishing for s > [16,21]. Because of the inversion 
symmetry of (1) [/(— x) = —/(#)], for each PO with characteristics (w p , A p ) there exists 
a PO with characteristics {—w p , A p ). This implies that d(~ 1 (P, s)/d(3\ j3=0 s>0 = 0. Using 
this fact and taking the limit of s — > + , we obtain the following expression for D: 

D = lim S^("')^a(~') (18) 



where 



The Hannay-Ozorio-de- Almeida uniformity sum rule [21,29] implies that lim ri ^ 00 (?(n) = 1. 
Then, the limit in (15) exists only if D WJ ^(n) converges to D. For sufficiently large n, D W \(n) 
in (17) is just the average of Wp/(2n) (p G U n ) weighted by the stability factor |A P | _1 . Such 
approximations to D have been used in previous works [15,18]. 

Formula (17) reduces to (4) in the case of uniform hyperbolicity. Consider, however, an 
equivalent expression for (17): 

avaw = 2n L n) Yj n m s m w2 i ( is ) 

w 

where N(n) and N w (n) are defined as in (4) and 
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S w (n) 



g(n)N w (n) 



nN{n) 



E 



(1-IApI- 1 ) 2 ' 



I A; 




(19) 



The quantity (19) is a natural restriction of (16) to the subset of POs with given winding 
number. If we now assume, in analogy to lim^oo g{n) = 1, the uniformity sum rules at 
fixed winding number w, 



formula (18) may reduce essentially to (4) also in cases of nonuniform hyperbolicity. We 
found much numerical evidence for the validity of (20) in our system. Part of this evidence is 
presented below. In general, the origin of (20) can be understood as follows. In our notation, 
one has the approximate relation for n 3> 1 (see, e.g., Appendix B in Ref. [30]): 



where r (an integer) is the repetition index. The left-hand side of (21) gives the probability 
distribution for a generic chaotic ensemble to diffuse a "distance" |w| in "time" n. Now, 
as n — > oo, there should be no essential difference between such an ensemble and the PO 
ensemble U n . The probability distribution above is then expected to be approximately equal 
to N w {n)/N{n) provided is not too close to the maximal value of \w p \ [10]. We use this 
in (21), keeping only the dominant terms (r = 1) on the right-hand side. Recalling also the 
definition (19) and the uniformity sum rule lim n ^ 00 g(n) = 1, Rel. (20) is obtained. 

In our numerical calculations we have used the relevant POs on the cylinder with peri- 
ods n < 14, computed as described in Sec. II. The following PO quantities were calculated 
accurately for k = 0.086A;, k = 1,..., 11: (a) S w (n) for all possible values of w > [note 
that S- W (n) = S w (n) from inversion symmetry], (b) The curvature-expansion (CE) approx- 
imation to D, DcE(n), obtained by using in (5) the expansion (7) truncated after the first n 
terms, (c) D^p^{n) [weighted-average formula (17)]. (d) D(n) [nonweighted-average formula 
(4)]. The PO results for D were compared with standard ones obtained from (2) by choosing 



lim S w {n) = 1, 



(20) 



n 




(21) 
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S as the entire unit torus. For this ensemble, which is obviously invariant under the torus 
map, one has the exact expansion [2] 



where Cj = (f(x )f(xj)) £ are the force-force correlations for (8). These correlations were 
calculated very accurately for j < 30 by a sophisticated integration of f(xo)f(xj) over the 
unit torus. In general, we found that D £ (m) in (22) converges rapidly to D due to the fast 
decay of Cj. For example, D £ (20) differs from both D £ (10) and D £ (30) by no more than 
0.05% for all the values of k considered. In what follows, D £ (20) will serve as our "standard" 
value D s for D. 

In Tables III-V we list the four quantities above for three representative values of k; the 
corresponding value of D s is also given. Table VI shows the quantities 5*0(14), Dce(14), 
D WA (14), -D(14), and D s for all the 11 values of k considered. The results in Tables III-V, 
as well as similar results for the other values of k, indicate that D CE (n), D w ^(n), and D(n) 
start to converge, in general, for n > 8. These quantities vanish for n < 3 since all the POs 
with period n < 3 have w p = 0, even for k = 0.946. Table VI shows a very good agreement 
between Z>wa(14) and Ds for all values of n. This agreement is generally better than that 
between _D CE (14) and D s . The relative difference between D WA (14) and D s ranges from 
0.04% to 0.4% while that between D CE (U) and D s ranges from 0.16% to 1.7%. 

Table VI also shows that the relative difference between the values of -D(14) and D$ for 
k < 0.43 is not larger than 2% despite the fact that the hyperbolicity for k < 0.43 can be 
considerably nonuniform, see Fig. 2. To understand this, consider the behavior of S w (n) 
for k < 0.43 in Tables III and IV. While the convergence of S w (n) to 1 for \w\ < 2 is quite 
evident, 53(14) and 5 4 (14) are not sufficiently converged due to the relatively small number 
of POs with w = 3, 4. Precisely because of this last fact, however, the effect of £±3(14) — 1 
and S±4(14) — 1 in (18) is not significant, leading to only a small difference between -D(14) 




m > 1, 



(22) 
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and D W a(14) (or D s ). As k is increased, a larger value of \w\ (\w\ = 5) appears for k > 0.516 
and the hyperbolicity becomes more nonuniform. Then, since the total number N(n) of POs 
is constant (does not depend on /t), the convergence of S w (n) is expected to deteriorate for 
all w. This is, in fact, confirmed by all our numerical data (see the worse case for k = 0.946 
in Table V) with the exception of w — 0. As Table VI shows, 5*0(14) remains remarkably 
well converged for all values of k. We shall attempt to find an explanation of this and other 
facts in a future work. 



IV. CONCLUSIONS 

In this paper, chaotic diffusion on periodic orbits (POs) was studied for the perturbed 
Arnol'd cat map on the cylinder, in a relatively large range of perturbations k corresponding 
to a "structural-stability" regime of the system on the torus. Numerical evidence indicates 
that this regime extends, at least approximately, significantly beyond the Anosov bound 
k c Ri 0.437, i.e., at least up to k ~ 1. This extension, which was already noticed in a 
quantum-chaos context for a different perturbed cat map [22], is further supported by the 
very good agreement between the PO and standard results for D also for k c < k < 0.946, 
where only POs topologically conjugate to the k = POs are used. In the absence of 
bifurcations, the variation of D with k is totally due to the change of the characteristics 
(w p , A p ) of a constant number of POs. Thus, the case studied in this paper is basically 
different from that considered in Ref. [17], i.e., standard maps in a strong-chaos limit. In 
the latter case, bifurcations of low-period (n — 1 and n = 2) POs are the main cause for 
the relevant variation of D with the parameter. The usual distinction between fundamental 
(n = 1) and curvature (n > 1) terms [19] is not felt in our case since w p = for n < 3 and 
for all the values of k considered. The convergence of the PO results to D generally starts 
only for order (or period) n > 8. 

As in the case of the Lorentz gas [18], the most accurate PO results are obtained by using 
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the weighted-average formula (17). In general, by expressing this formula in the winding- 
number representation (18), it becomes clear that the effect of a nonuniform hyperbolicity is 
completely captured by the basic quantities (19). Since these quantities satisfy the unifor- 
mity sum rules (20), the manifestation of this effect is essentially restricted, for sufficiently 
large n, to the "tail" of the distribution (21) (\w\ large relatively to the maximal value of 
\w p \). Here the discrepancy between N w (n)/N(n) and (21) leads to a value of S w (n) which 
is not well converged. We have shown that the effect of nonconverged values of S w (n) may 
be insignificant also in cases of a considerably nonuniform hyperbolicity. Then, formula (17) 
reduces essentially to the nonweighted-average formula (4). 
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Table I. ( 1 (0, 0) to order n — 1, 14 of curvature expansion for several values of re. 



n 


C-^o, o) 




re = 0.086 


re = 0.258 


re = 0.43 


re = 0.602 


re = 0.774 


re = 0.946 


1 


-0.36790 


-0.34306 


-0.32172 


-0.30313 


-0.28676 


-0.27220 


2 


-0.29773 


-0.31199 


-0.33060 


-0.35530 


-0.38853 


-0.43399 


3 


-0.17097 


-0.17537 


-0.17565 


-0.17293 


-0.16915 


-0.16734 


4 


-0.08734 


-0.09026 


-0.09120 


-0.08990 


-0.08641 


-0.08143 


5 


-0.04179 


-0.04339 


-0.04390 


-0.04242 


-0.03686 


-0.02308 


6 


-0.01918 


-0.02003 


-0.02043 


-0.01995 


-0.01786 


-0.01328 


7 


-0.00856 


-0.00899 


-0.00925 


-0.00909 


-0.00795 


-0.00478 


8 


-0.00374 


-0.00395 


-0.00410 


-0.00408 


-0.00360 


-0.00222 


9 


-0.00342 


-0.00688 


-0.00995 


-0.01261 


-0.01476 


-0.01621 


10 


-0.00002 


0.00104 


0.00185 


0.00249 


0.00306 


0.00375 


11 


0.00025 


0.00130 


0.00237 


0.00350 


0.00480 


0.00643 


12 


-0.00002 


0.00020 


0.00042 


0.00059 


0.00073 


0.00087 


13 


0.00019 


0.00061 


0.00097 


0.00125 


0.00147 


0.00164 


14 


0.00012 


0.00038 


0.00062 


0.00083 


0.00099 


0.00103 



Table II. Number N(n) of relevant POs of period n, having no point oni = —0.5. 



n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


N(n) 


1 


2 


2 


10 


24 


48 


120 


270 


568 


1500 


3600 


8543 


20880 


50700 
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Table III. S w (n), D CE (n), D WA (n), and D(n) for « = 0.258; D s = 0.04865. 



n 


S (n) 


Si(n) 


S 2 (n) 


S 3 (n) 


S 4 (n) 


-Dce(^) 


D WA (n) 


D{n) 


1 


1.52221 










0.00000 


0.00000 


0.00000 


2 


1.19039 










0.00000 


0.00000 


0.00000 


3 


1.05180 










0.00000 


0.00000 


0.00000 


4 


1.01987 


1.02505 








0.02724 


0.05125 


0.05000 


5 


1.03504 


0.97134 








0.03261 


0.04047 


0.04167 


6 


1.02088 


0.98225 








0.03350 


0.03752 


0.03819 


7 


1.01632 


0.96115 


1.13505 






0.04102 


0.05023 


0.04881 


8 


1.02287 


0.96467 


1.06022 






0.04410 


0.04795 


0.04768 


9 


0.99896 


0.99316 


1.04018 






0.04487 


0.04762 


0.04695 


10 


1.01179 


0.98401 


1.01905 


1.24670 




0.04648 


0.04807 


0.04760 


11 


1.00966 


0.98787 


1.00285 


1.17419 




0.04726 


0.04862 


0.04811 


12 


1.00935 


0.98897 


1.00102 


1.11421 




0.04751 


0.04819 


0.04776 


13 


1.00981 


0.98864 


0.99965 


1.08830 


1.40865 


0.04801 


0.04873 


0.04826 


14 


1.00658 


0.99180 


0.99950 


1.05561 


1.30305 


0.04824 


0.04861 


0.04815 
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Table IV. S w (n), D CE (n), D WA (n), and D(n) for k = 0.43; D s = 0.05380. 



n 


S (n) 


Si(n) 


S 2 (n) 


S 3 (n) 


S4 




-Dce(^) 


D WA (n) 


D{n) 


1 


1.47432 












0.00000 


0.00000 


0.00000 


2 


1.21592 












0.00000 


0.00000 


0.00000 


3 


1.04705 












0.00000 


0.00000 


0.00000 


4 


1.03876 


0.99752 










0.02630 


0.04988 


0.05000 


5 


1.06682 


0.92808 










0.03128 


0.03867 


0.04167 


6 


1.02173 


0.95875 


1.27251 








0.03923 


0.05429 


0.05208 


7 


1.02423 


0.95049 


1.13318 








0.04692 


0.05640 


0.05476 


8 


1.02965 


0.95429 


1.07500 








0.04871 


0.05239 


0.05185 


9 


0.99527 


0.98363 


1.06576 


1.74814 






0.05076 


0.05433 


0.05184 


10 


1.02227 


0.97717 


0.98799 


1.35259 






0.05216 


0.05339 


0.05253 


11 


1.01290 


0.98137 


1.00161 


1.23369 






0.05265 


0.05426 


0.05308 


12 


1.01111 


0.98476 


1.00131 


1.13098 


2.06872 


0.05294 


0.05318 


0.05216 


13 


1.00872 


0.98746 


0.99941 


1.08594 


1.59230 


0.05323 


0.05372 


0.05271 


14 


1.00869 


0.98753 


0.99789 


1.07197 


1.43? 


^66 


0.05338 


0.05378 


0.05273 
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Table V. S w (n), D CE (n), D WA (n), and D(n) for k = 0.946; D s = 0.07200. 



n 


S (n) 


Si(n) 


S 2 (n) 


S 3 (n) 


S 4 (n) 


S 5 (n) 


-Dce(^) 


D WA (n) 


D{n) 


1 


1.37400 














0.00000 


0.00000 


0.00000 


2 


1.42345 














0.00000 


0.00000 


0.00000 


3 


1.03401 














0.00000 


0.00000 


0.00000 


4 


1.19364 


0.77294 












0.01918 


0.03865 


0.05000 


5 


1.05191 


0.75048 


2.07839 










0.05413 


0.10055 


0.07500 


6 


1.14789 


0.70455 


1.63690 










0.06132 


0.06993 


0.06250 


7 


1.06410 


0.72943 


1.60419 










0.06600 


0.08195 


0.06667 


8 


1.12569 


0.80405 


1.16902 


1 


83040 






0.06826 


0.06781 


0.06389 


9 


0.94494 


0.90200 


1.33077 


1 


93640 






0.07233 


0.08259 


0.06494 


10 


1.00057 


0.98694 


0.91645 


1 


64065 






0.07232 


0.06972 


0.06433 


11 


0.96412 


0.96616 


1.09809 


1 


32642 


1.94196 




0.07138 


0.07231 


0.06485 


12 


1.01872 


0.94809 


1.02208 


1 


19621 


2.92740 




0.07175 


0.07120 


0.06448 


13 


0.97018 


0.95294 


1.10901 


1 


09196 


2.04263 


3.55700 


0.07181 


0.07262 


0.06507 


14 


1.00396 


0.94759 


1.04419 


1 


11549 


1.90135 


2.53191 


0.07188 


0.07180 


0.06505 
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VI. 50(14), £>ce(14), £>wa(14), D(U), and D s for all the 11 values of 



K 


50(14) 


£>ce(14) 


Ava(14) 


D(U) 


£s 


0.086 


1.00193 


0.04350 


0.04377 


0.04356 


0.04388 


0.172 


1.00431 


0.04551 


0.04618 


0.04579 


0.04622 


0.258 


1.00658 


0.04824 


0.04861 


0.04815 


0.04865 


0.344 


1.00752 


0.05070 


0.05109 


0.05046 


0.05116 


0.43 


1.00869 


0.05338 


0.05378 


0.05273 


0.05380 


0.516 


1.00807 


0.05589 


0.05642 


0.05506 


0.05656 


0.602 


1.00058 


0.05841 


0.05938 


0.05736 


0.05942 


0.688 


0.99577 


0.06220 


0.06228 


0.05921 


0.06240 


0.774 


1.00191 


0.06531 


0.06548 


0.06102 


0.06551 


0.86 


1.00014 


0.06838 


0.06848 


0.06297 


0.06873 


0.946 


1.00396 


0.07188 


0.07180 


0.06505 


0.07200 
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FIGURES 



Fig. 1. Mixed phase space for the perturbed Arnol'd cat map [(1) with (8)] for k — 2. 

Fig. 2. Distributions of the Lyapunov exponents \ p = ln(|A p | _1 )/n for all the primitive 
POs of period n = 14 of the perturbed Arnol'd cat map for k = 0.258 (filled circles), k = 0.43 
(crosses), and k = 0.946 (triangles). Each distribution was calculated by dividing the full 
range of A p into 100 equal intervals and counting the number AiV of values of \ p in each 
interval. 
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